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Abstract 

We study the polarization observables of J/ip hadroproduction at next-to-leading 
order within the factorization formalism of nonrelativistic quantum chromodynam- 
ics. We complete the present knowledge of the relativistic corrections by also provid- 
ing the contribution due to the intermediate 3 -fj 8 ' color-octet states at this order, 
which turns out to be quite significant. Exploiting the color-octet long-distance 
matrix elements previously extracted through a global fit to experimental data of 
unpolarized J/ip production, we provide theoretical predictions in the helicity and 
Collins-Soper frames and compare them with data taken by CDF at Fermilab Teva- 
tron I and II and by ALICE at CERN LHC. The notorious CDF J/tp polarization 
anomaly familiar from leading-order analyses persists at the quantum level, while 
the situation looks promising for the LHC, which is bound to bring final clarification. 
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While the overly successful experiments at the LHC are closing in on the Higgs boson 
and are systematically searching for signals of physics beyond the standard model (SM), 
we must not be carried away losing track of a longstanding, unresolved puzzle in quantum 
chromo dynamics (QCD), the otherwise well-established SU(3) gauge theory of the strong 
interactions, right in the core of the SM. In fact, despite concerted experimental and 
theoretical efforts ever since the discovery of the J/ip meson in the November revolution of 
1974 (The Nobel Prize in Physics 1976), the genuine mechanism underlying the production 
and decay of heavy quarkonia, which are QCD bound states of a heavy quark Q = c,b 
and its antiparticle Q, has remained mysterious. 

Nonrelativistic QCD (NRQCD) (Tj endowed with an appropriate factorization theo- 
rem, which was conjectured in a seminal work by Bodwin, Braaten, and Lepage [2] and 
explicitly proven through next-to- next-to-leading order [3], arguably constitutes the most 
probable candidate theory at the present time. This implies a separation of process- 
dependent short- distance coefficients, to be calculated perturbatively as expansions in 
the strong-coupling constant a s , from supposedly universal long-distance matrix elements 
(LDMEs), to be extracted from experiment. The relative importance of the latter can be 
estimated by means of velocity scaling rules, which predict each of the LDMEs to scale 
with a definite power of the heavy-quark velocity v in the limit v < 1. In this way, the 
theoretical predictions are organized as double expansions in a s and v. A crucial feature 
of this formalism is that the QQ pair can at short distances be produced in any Fock state 
n = 2S + 1 L / y with definite spin S, orbital angular momentum L, total angular momentum 
J, and color multiplicity a = 1,8. In particular, this formalism predicts the existence of 
intermediate color-octet (CO) states in nature, which subsequently evolve into physical, 
color-singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the limit 
v — > 0, the traditional CS model (CSM) is recovered in the case of S-w&ve quarkonia. In 
the case of J/ip production, the CSM prediction is based just on the 3 Si CS state, while 
the leading relativistic corrections, of relative order 0(v 4 ), are built up by the x 4 8] , 3 Sf ] , 
and 3 Pj 8] (J = 0, 1, 2) CO states. 

The CSM is not a complete theory, as may be understood by noticing that the NLO 
treatment of P-wave quarkonia is plagued by uncanceled infrared singularities, which are, 
however, properly removed in NRQCD. In a way, NRQCD factorization [2], appropriately 
improved at large transverse momenta by systematic expansion in powers of tt?q/p\ 
[3], is the only game in town, which makes its experimental verification such a matter of 
paramount importance and general interest [5]. 

The experimental test of NRQCD factorization [2] has been among the most urgent 
tasks on the agenda of the international quarkonium community for the past fifteen years 
[5] and, with high-quality data being so copiously harvested at the LHC, is now more 
tantalizing than ever. The present status of testing NRQCD factorization in charmonium 
production is as follows. As for the unpolarized J/ip yield, NRQCD factorization has 
recently been consolidated at NLO by a global fit to the world's data of hadroproduction, 
photoproduction, two-photon scattering, and e + e~ annihilation [6], which successfully 
pinned down the three CO LDMEs in compliance with the velocity scaling rules and 
impressively supported their universality (see Table 1 in Ref. [7]). On the other hand, the 
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Tevatron II |8J data alone can just fix two linear combinations of the three CO LDMEs 
[9][T0] . and the fit results of Ref. [9] are incompatible with Ref. [6], as discussed in Ref. [TT] . 
As for the J/ip polarization observables, a complete NLO NRQCD analysis has so far 
only been performed for photoproduction [7j. The agreement with recent data from 
DESY HERA [12] was found to be satisfactory at sufficiently large values of pt, where 
nonperturbative soft-gluon effects are negligible, provided that diffractive events were 
eliminated by an appropriate acceptance cut. The overall \ 2 achieved at NLO in NRQCD 
turned out to be more than a factor of 2 below the value obtained at NLO in the CSM 
[T3] . However, the case for NRQCD factorization appeared here to be not as convincing 
as for the unpolarized J /if) yield [6]. 

In this Letter, we take the next — and possibly decisive — step in the worldwide endeavor 
to put NRQCD factorization to the test and scrutinize polarized J /if) hadroproduction at 
NLO in NRQCD. Measurements by CDF in pp collisions with cm. energy yfs = 1.8 TeV 
[2] and y/s = 1.96 TeV [TB] at the Tevatron and by ALICE [TB] in pp collisions with 
a/s = 7 TeV at the LHC are waiting to be rigorously interpreted. Previous NLO analyses 
were confined to the CSM [17] or included only the ^q 8 ' and 3 S® contributions [18], which 
are deducible by well-established techniques. We close this gap by providing also the 3 Pj 8 ^ 
contributions at NLO, which are expected to be significant numerically. Their calcula- 
tion is far more intricate because the applications of the respective projection operators 
to the short-distance scattering amplitudes produce particularly lengthy expressions in- 
volving complicated tensor loop integrals and exhibiting an entangled pattern of infrared 
singularities. Technical details will be presented in a forthcoming publication. 

The polarization of the J /if) meson is conveniently analyzed experimentally by mea- 
suring the angular distribution of its leptonic decays, which is customarily parametrized 
using the three polarization observables Xg, A^, and Xg^, as 

W(9, 0) oc 1 + Xg cos 2 6 + A sin 2 6 cos(20) + \ e<t> sin(20) cos 0, (1) 

where 9 and are respectively the polar the azimuthal angles of / + in the J /if) rest 
frame. This definition depends on the choice of coordinate frame. In the experimental 
analyses [T4|15|16] . the helicity and Collins-Soper frames were employed, in which the 
polar axes point in the directions of — (p p + pp) and p p /\p p \ — Pp/\pp\, respectively. The 
values Xg = 0, +1, — 1 correspond to unpolarized, fully transversely polarized, and fully 
longitudinally polarized J /if) mesons, respectively. In Refs. [Ll|T5] . Xg is called a. Working 
in the spin density matrix formalism and denoting the z component of S by i, j = 0, ±1, 
we have 

_ dan - da 00 \/2Redcr 10 

Ae ~ 1 TTt ' ^ — 1 TTt — ' Ke <t> ~ i — n ' 

d(j\\ + ao"oo ao"ii + cLctqq aan + acroo 

where da^ is the ij component of the pp — > J /if) + X differential cross section. An 
expression of da^ in terms of parton density functions (PDFs) and partonic spin density 
matrix elements may be found in Eq. (3) of Ref. [7J. 

In our numerical analysis, we adopt the CO LDME values from Table I of Ref. [7J along 
with (C J ^( 3 ,S1 1] )) = 1.32 GeV 3 [19J, take the charm-quark mass, which we renormalize 
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Figure 1: (color online) NLO NRQCD predictions (solid lines) for X e and A^ as func- 
tions of pt in the helicity and Collins-Soper frames including theoretical uncertainties 
(shaded/yellow bands) compared to CDF [T4"yT5] and ALICE [16] data. For compari- 
son, also the NLO CSM (dot-dashed lines) predictions including theoretical uncertainties 
(hatched/blue bands) as well as the LO NRQCD (dashed lines) and LO CSM (dotted 
lines) ones are shown. 

according to the on-shell scheme, to be m c = 1.5 GeV, and use the one-loop (two-loop) 
formula for ai n/ ^(/i r ), with rif = 4 active quark flavors, at LO (NLO). As for the proton 
PDFs, we use the CTEQ6L1 (CTEQ6M) set [20J at LO (NLO), which comes with an 
asymptotic scale parameter of Aqq D = 215 MeV (326 MeV). Our default choices for the 
MS renormalization, factorization, and NRQCD scales are fi r = [if = rriT and fi\ = 

m c , respectively, where my = \Jp T + 4m^ is the J/ip transverse mass. The theoretical 
uncertainty due to the lack of knowledge of corrections beyond NLO is estimated by 
varying /i r , [if, and fi\ by a factor 2 up and down relative to their default values. In 
our NLO NRQCD predictions, we must also include the errors in the CO LDMEs, which 
reflect the errors on the experimental data included in the fit. To this end, we determine 
the maximum upward and downward shifts generated by independently varying their 
values according to Table I in Ref. [7] and add the resulting half-errors in quadrature to 
those due to scale variations. 

In Fig. [TJ we confront our NLO NRQCD predictions for Xg and as functions of 
Pt in the helicity and Collins-Soper frames with the measurements by CDF [14115] and 
ALICE [16]. Since the cross section ratios in Eq. (T2]) are very insensitive to the precise 
value of y/s, we may safely overlay the data from y/s = 1.8 TeV [13] with the predictions 
for y/s = 1.96 TeV. For comparison, also the LO NRQCD as well as the LO and NLO 
CSM predictions are shown. In order to visualize the size of the NLO corrections to the 
hard-scattering cross sections, the LO predictions are evaluated with the same LDMEs. 
As in Ref. [6] , we do not consider the range pt < 3 GeV, where nonperturbative soft-gluon 
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radiation invalidates a fixed-order treatment. We observe that, in all the cases considered, 
the inclusion of the NLO corrections has a considerably less dramatic effect in NRQCD 
than in the CSM, where the normalizations and shapes of the various distributions are 
radically modified. This indicates that the perturbative expansion in a s converges more 
rapidly in NRQCD than in the CSM. A similar observation has already been made for 
photoproduction [7J. Comparing the LO NRQCD prediction for Xg in Fig. [Qa) with 
the one referring to direct J/ip production in Fig. 1(b) of Ref. [21], we encounter marked 
differences in normalization and shape, which mainly reflect the progress in our knowledge 
of the CO LDMEs during the dozen years between Ref. [21] and [6]. Comparing Figs. [T^a) 
andH^b), we observe that, for px ^> 2m c , the results are fairly stable when passing from 
the Tevatron to the LHC and from the central (\y\ < 0.6) to the forward (2.5 < y < 4) 
rapidity region, as long as one stays in the helicity frame. However, switching from the 
helicity to the Collins-Soper frame has a radical impact on the various px distributions, 
as the comparison of Figs. QJb) and []Jc) reveals. The most striking effect appears for 
Xg, where the NLO CSM and NRQCD results are, roughly speaking, inverted. As for 
A^, the Collins-Soper frame clearly outperforms the helicity frame with regard to the 
power of NRQCD versus CSM discrimination. As expected, the theoretical uncertainties 
due to scale variations steadily decrease as the value of px increases, which just reflects 
asymptotic freedom. 

Let us now compare experiment and theory. The Tevatron I data for \q [H] systemat- 
ically fall below the NLO NRQCD prediction and, unlike the latter, exhibit a downward 
trend for p^^lO GeV, but their errors are too large for a firm conclusion. This is quite 
different for the Tevatron II data [15] , which are rather precise and indicate that the J / ip 
mesons are essentially unpolarized in the helicity frame, while NRQCD clearly predicts 
transverse polarization, with purity and precision steadily increasing with the value of px- 
A few caveats are in order here. While the CDF data cover prompt production, including 
the feed-down from the heavier Xcj and tp' mesons, our predictions refer to direct produc- 
tion, devoid of feed-down. However, the inclusion of feed-down was found to have a minor 
effect on the px distribution of Xg at LO in NRQCD [21] . which is expected to carry over 
to NLO. For px < 12 GeV, the measurements from Tevatron runs I [2] and II [TS] are 
mutually incompatible, a feature that has never been satisfactorily clarified by the CDF 
Collaboration. Fortunately, the four LHC experiments are in a position to measure the 
J/ip polarization with unprecedented precision, and ALICE has already presented first 
results for Xg and A^, [16], both in the helicity and Collins-Soper frames. We compare 
them with our NLO NRQCD predictions in Figs. [T^b) and[]Jc), respectively, leaving out 
one data point in the px bin 2-3 GeV in each of the four cases. The degree of agreement 
is encouraging. All data points are at most one standard deviation away from the the- 
oretical error bands; four of them even overlap with the latter. The agreement tends to 
improve with increasing value of px, as the influence of uncontrolled corrections beyond 
NLO fades. Except in the \ow-px bin of Xg in the helicity frame, the data clearly favor 
NRQCD over the CSM at NLO. However, also these comparisons have to be taken with a 
grain of salt, since, unlike our predictions, the ALICE data [IB] also include J/ip mesons 
from feed-down and B- meson decays. The fraction of J/ip mesons from 6-quark origin 
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Figure 2: (color online) da 00 /dpT, d<Tu/dp T , and dai_i/dp T for pp — > J/tj) + X with 
yfs = 7 TeV and 2.5 < ?/ < 4 in the helicity and Collins-Soper frames at NLO in NRQCD 
(solid lines) and their decompositions in the 3 Sf^ (dot-dashed lines), ^q 8 ' (dotted lines), 
3 Si (coarsely dotted lines), and 3 Pj 8 ^ (dashed lines) contributions. Negative values are 
marked red. 

was measured by the LHCb Collaboration for unpolarized production to be below 15% in 
the pt range considered here [22] , so that their omission should have an effect negligible 
against the theoretical uncertainty. 

In order to assess the relative importance of the individual cc Fock states n for Xg and 
in the helicity and Collins-Soper frames, we detail their contributions to daoo, dan, and 
dai-i in Fig. |2j Note that the unpolarized cross section is recovered as daQQ + 2dau, while 
do\-\ is an auxiliary quantity, which is entitled to take on either sign and receives no 
contribution from the ^q 8 ' channel. As anticipated above, the previously unknown 3 Pj 8 ^ 
contributions play a dominant role in this game, and their omission is bound to yield 
meaningless results. As is familiar from the unpolarized case [B], the 3 Pj 8 ' short-distance 
cross sections receive sizable NLO corrections that may even turn them negative beyond 
some values of pt, as may be inferred from Figs. [2]^b), (d), and (e) by recalling that our 
NLO value of (O j ^( 3 P^)) [6J is also negative. This is, however, unproblematic because 
a particular CO contribution represents an unphysical quantity depending on the choices 
of renormalization scheme and scale fi\ and is entitled to become negative as long as the 
full cross section remains positive. 
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At this point, we compare our results with the theoretical literature. We agreed with 
the LO NRQCD results given in Eqs. (B26)-(B62) of Ref. [23]. We are able to nicely 
reproduce the LO and NLO 3 S{ 1] and 3 Sf ] contributions to Xg shown as functions of pr in 
Fig. 3 (9) of the first (second) paper in Ref. [17] and Fig. 4 of Ref. [T8] . 

In conclusion, we presented the first complete NLO analysis of polarized J /if) hadropro- 
duction within NRQCD and, exploiting the knowledge of the CO LDMEs from a global 
fit of unpolarized J /if) production data [6], thus managed to establish certainty about 
that the latest CDF measurement [15] is in severe conflict with NRQCD factorization [2] . 
As emphasized in Ref. [24j, the measurement of Xe in the helicity frame alone provides 
only restricted access to the rich phenomenology of J/if> polarization and might become 
ambiguous in want of statistics, in particular if the patches of solid angle in which sim- 
ulations substitute data exceed a critical size. In this sense, the ALICE Collaboration 
[T6] has made a first step into a bright future of heavy-quarkonium physics at the LHC, 
allowing for ultimate tests of NRQCD factorization. 

This work was supported in part by BMBF Grant No. 05H09GUE and HGF Grant 
No. HA 101. 

Note added. After the submission of this Letter, a preprint |25j appeared, in which 
Xe in the helicity frame is studied at NLO in NRQCD for J /if) hadroproduction at the 
Tevatron. Adopting the LDME values specified in the first line of Table I in Ref. [25], we 
are able to reproduce the numerical results of Ref. [25J. Unfortunately, including the CDF 
data of Ref. [14115] on top of those from Ref. [8], already used in Ref. [9], still does not 
allow for an independent determination of the three CO LDMEs, and the result is violently 
incompatible with the global fit of Ref. [6]. In particular, this set of CO LDMEs, as well 
as the two alternative ones specified in that table, are found [11] to yield NLO-NRQCD 
predictions for photoproduction which fail to describe the HERA data. As a consequence, 
the conclusions of Ref. [25J with regard to the compliance of the CDF measurements of 
Xe in the helicity frame [T4"yT5] with NRQCD factorization at NLO are quite opposite to 
ours. 
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